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SUMMARY 


A study of the energy absorption mechanisms and fluid dynamic 
considerations for efficient conversion of high power laser radiation into a 
high velocity flow is presented. Although the contents of this study are 
applicable to the laser propulsion concept, they are not unique to this 
application and are useful for analysis of any system which requires the 
efficient conversion of laser energy to kinetic energy. The objectives are 
(1) to determine the most effective absorption mechanisms for converting 
laser radiation into translational energy, and (2) examine the require- 
ments for transfer of the absorbed energy into a steady flow which is stable 
to disturbances in the absorption zone. The first part of the study consists 
of a review of inverse Bremsstraklung, molecular and particulate absorp- 
tion mechanisms. 

Inverse Bremsstrahbing offers the advantage of spectral continuum 
absorption along with the rapid conversion of absorbed energy into trans- 
lational modes. However, some mechanisms must be applied to provide an 
initial electron concentration, and until the gas becomes fully ionized, the 
absorption coefficient is strongly temperature dependent. The latter may 
induce unstable heating. A review of experimental and theoretical work in 
electron- neutral and electron- ion Bremsstrahlung is presented. 

Molectilar absorption offers the advantage of significant absorption 
at low temperatures plus the possible weak dependence of the absorption co- 
efficient upon temperature. In general, however, molecular absorption 
requires coincidence between wavelengths of lasing and absorbing transitions, 
may ’*bleach” at large laser fluxes, requires additional techniques for con- 
verting vibrational energy to translational energy and may lose its efficiency 
for absorption via chemical reactions. It is proposed that the first and 
second of the above disadvantages of molecular absorption can be mitigated 
by high pressures, and high pressure molecular absorption is modeled. 

In addition, an analysis of a potential system for converting the absorbed 
energy into translation by using atomic species to vibrationally deactivate 
the absorbing molecule is presented. 

Particulates offer the advantages of continuum absorption, rapid 
conversion of laser energy to translational energy and high conversion effi- 
ciency since only a small percentage of total absorbed energy remains in 
the particle. Modeling of the gas heating characteristics suggests that 
optimization occurs for very small particles at mass loadings of the order 
of lO""^ gms/cm^. As the particles approach. the boiling point other mech- 
anisms such as molecular absorption and Bremsstrahlung can become the 
dominant absorption mechanisms, and v^aporization and radiative losses from 
the particles may decrease their absorption efficiency. 
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The second part of the study consists of steady flow and stability 
considerations for conversion of the laser power to a high velocity flow in 
a nozzle configuration. The quasi-one -dimensional flow through a nozzle 
is formxilated under the assumptions of perfect gas, instantaneous conver- 
sion of absorbed laser energy to temperature (equilibrium flow) and an 
absorption coefficient proportional to density and temperature raised to 
arbitrary powers* For a specified nozzle configuration, predictions of Mach 
number, temperature, density and exhaust velocity are presented as a 
function of optical depth and the ratio of the las ir power to inlet flow power. 
A "local" stability analysis is performed and a typical stability map of 
disturbance wavenumber versus nozzle position is presented. Neutral 
stability contours provide an indicator for proper design of nozzles with 
stable absorption zones. 
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1. INTRODUCTION 

The success of any system designed to convert high power laser 
radiation into directed kinetic energy depends upon the ability of the system 
to perform this conversion efficiently. Neither the existing research on the 
use of remotely stationed lasers to power a rocket engine^ nor the work 
of Hertzberg et al^ on the thermodynamics of a photon engine have a priori 
addressed the fundamental issue of the most efficient process by which a 
flowing gas can absorb laser radiation and convert it to useful work. The 
research in laser propulsion^# 2 relied upon the creation of electrons via 
breakdown to initiate absorption by inverse Bremsstrahlung and subsequent 
conversion of the heated gas into kinetic energy to obtain high specific 
impulse. However, the experiments^ illustrated a posteriori that (for solid 
propellants) without some tailoring of the propellant composition, the con- 
version process yields the initiation of absorption waves in the gas. 2 The 
result was unsteady flow behavior and low power conversion efficiency. 

No attempts were made to determine whether other absorbers (e.g. molecular 
or particulate) would provide a stable and efficient method for converting 
laser radiation into directed kinetic energy. 

What is necessary for any application which requires the conversion 
of coherent laser radiation into flow work are answers to the following 
questions: 

(1) What are the best absorbers to make this conversion from 
radiation to translational energy most efficient? 

(2) Is the absorption of laser radiation by a flowing gas a 
steady process where the gas enters the absorbing region 
in one state and exits in another? In addition, upon 
leaving the absorption zone, has the absorbed energy 
been stored in an excited molecular state to be converted 
to translational energy downstream? 

(3) Is the conversion process stable to disturbances in the 
absorbing region? 

This report is aimed at providing the first step toward answering 
the above questions. The report is essentially divided into two major parts. 
The next three sections art* devoteci to a review of absorption mechanisms to 
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convert laser energy to translational energy. Section 2 consists of a 
review of inverse Bremsstrahlung absorption. Section 3 discusses the 
properties of molecular absorption and presents an analysis of a proposed 
molecular absorption system at high pressure. Section 4 is devoted to 
a discussion of radiative absorption via particulates. The last section 
describes the second part of our program and presents fluid dynamic con- 
siderations for converting the absorbed laser energy to a high velocity 
flow. Steady solutions for the quasi-one-dimensional nozzle flow of a 
coherent radiation absorbing gas are presented. In addition, a stability 
analysis is described and applied to the nozzle flow to determine under 
what conditions the conversion process is stable to disturbances in the 
absorbing region. 
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Z. INVERSE BREMSSTRAHLUNG 
2. 1 INTRODUCTION 

When a free electron is in a close encounter with an atom or 
molecule it will undergo an acceleration caused by the potential field betueen 
the collision pair. In the case of collision with an ion the dynamics of the 
collision is dominated by the long range coulombic (or screened coulombic) 
field whereas in the collision with a neutral species the resulting induced 
fields^ dipole or higher, are short range in nature. It is well known from 
classical electrodynamics that an acv elerating electron radiates continuously, 
and in the process loses part of its kinetic energy, thus, slowing down. For 
this reason the resulting radiation is called Bremsstrahlung (braking radia- 
tion) and the reverse absorption process is termed inverse Bremsstrahlung. 
The radiation resulting from electron- ion collisions is often called 
Kramers* radiation while that resulting from electron- neutral collisions 
is termed neutral Bremsotralilung. 

Bremmstrahlung radiation is of interest in many areas including 
astrophysical applications, laser heating of plasmas, and in the prediction 
of the infrared radiation from heated gases and has received considerable 
theoretical study. The interest of this work is in the prediction of the in- 
frared inverse Bremsstrahlung absorption coefficient for light gases, such 
as Hydrogen or Helium, at temperatures less than 10,000^K. Recent^* ^ 
detailed predictions and measurements of the absorption coefficients for these 
gases have been performed at temperatures above 10,000^K. In this instance 
the comparison between theory^ and experiment*^ was excellent. Un- 
fortunately, at lower temperatures, where the neutral Bremsstrahlung 
process dominates, there do not seem to be any experimental measure- 
ments for H and He and theoretical estimates must be relied upon. 

Although a definitive review is not justified under the present effort, 
a general discussion of the theoretical work on Kramers radiation and neutral 
Bremsstrahlung will be given below. Perhaps the most detailed discussion 
of the semi-classical treatment of Kramers radiation and the related free 
bound processes may be found in the Handuch d. Physik article by Finkelnburg 
and Peters.® Both Biberman and Norman^ and Johnston^® have presented 
surveys of the theoretical work performed in this area up to 1967. A number 
of more recent efforts are referenced in Ref, 5 and 6. The available 
measurements of neutral Bremsstrahlung cross-sections at infrared wave- 
lengths are limited and to our knowledge have never been reviewed. 
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In this section the data base will be compared with theoretical predictions 
wherever possible. 

2. 2 KRAMERS RADIATION (ELECTRON-ION BREMSSTRAHLUNG ) 

A free electron by itself cannot absorb a photon and conserve both 
momentum and energy, however, if a third body is available to pick up the 
recoil momentum, the process is allowed. Kramers was the first to 
develop a semi- classical model for this phenomenon for the case of an 
electron interacting with a hydrogenic ion in a pure coulomb field. His 
result for the absorption coefficient was 

+ hv/kT j 

\ n. (e - 1), cm' (2. 1) 


where is the electron-ion absorption cross-section, n^ and n^, the 
electron and ion densities respectively, and v the frequency of the laser. 

The bracketed term represents a product of two distinct effects. The first 
being a factor cf (l-e“^^' ) which is a correction term for stimulated ennis- 

sion, and the second, a factor of which includes the effect of bou£id-free 

absorption from high lying electronic states. This latter term was derived^^ under 
the assumption that the electron and electronic temperatures are in equilibrium 
and should not be used indiscriminately. In anv event this latter term is un- 
important for hv/kT« 1. 

The quantity 0^^ is given by 


4 / 2tt \ 5 

3 . 3 ' 

\ ' h c m V 

3. 7 X 10^ 


( 2 . 2 ) 


where Z is the effective nuclear charge, e is the unit charge, m is the 
electron mass, k is Boltzmanns constant and h is Planks constant. The 
basic assumptions in this derivation are a pure coulomb potential field 
(hydrogenic ions) and an electron velocity distribution defined by a 
Maxwellian distribution at the temperature T. 
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The quantum mechanical derivation was first considered by Gaunt. 

The results are complex but may be evaluated accurately for the case of a 
hydrogenic ion. The quantum mechanical solution is generally written in 
terms of the classical solution, Eq. (2.2), multiplied by a correction factor 
g (^* T), called the Gaunt Factor. There is no simple formulation for the 
Gaunt Factor, however, it has been tabulated extensively by Karzas and 
Latter^^ among others. For our conditions of interest, 0. 1 < hv 1 ev, 
1000 < T < 10,000°K, the Gaunt Factor is ~ 1. 3 ± 0.2. 

The quantum mechanical prediction for ions other than hydrogenic 
is considerably more complicated because of screening effects as well as 
distortion of the electronic cloud. However, the long range coulomb field 
can dominate the momentum absorption process for low energy photons 
resulting in cross sections similar to that for hydrogenic ions. This point 
is discussed by Johnston, and small first order corrections for species 
such as He have been calcvilated by Peach. 

2. 3 NEUTRAL BREMSSTRAHLUNG 

2. 3. 1 Theory 

Free-free absorption in the presence of a neutral atom will dominate 
Kramers radiation when the fractional concentration of ions is % 10“^. The 
basic complication in evaluating the cross sections for this process is in 
defining the interaction potential between the electron and the neutral species. 
This is particularly complicated in the case of an n electron atom or 
molecule and generally an effective spherically symmetric potential is 
assumed. 

Following Geltman^^ the differential cross section for the absorption 
of a photon of energy hv by an electron of initial wave vector k j resulting 
in the final wave vector kf is given by^^ 


2 

do e mvk 

en 

dn 


f , , /N , 

2 |< f |‘ • r(i >1 


(2.3) 


f 'K c 

where is the electric field polarization unit vector and | i> and | f > are 
Lmu 

states is 


the contimuum states partial wave expansion of the continuxim 

2-f ' (2AM).- — 


I 0 
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where iijj is the scattering phase shift and (k« r) is the radial wave 
function defined 




2 , 2 

d /dr - 


{»■ ^ 1 ) 


- U (r) + k 


In- 


(2.5) 


where U (r) is the potential function defining the electron- neutral interaction. 
The overall absorption cross section, is determined by integrating 

Eq. (2. 3) over all incident and final electron directions and then performing 
a Maxwell average. The resulting absorption coefficient is given by^^ 


hv/kT 


= Q 

en '^en 




"n<* 


), cm 


-1 


( 2 . 6 ) 


where n^ and are the electron and neutral densities respectively and 
the bracketed term is a correction for stimulated emission. 

Early evaluations ^ of exprossion(2.3) were carried out to first 
order, s wave only, and the resulting terms could be related to the elastic 
scattering or momentum transfer cross sections, thus eliminating the need 
for specification of the interaction potential. More recent calculations 
(refs. 16, 19, 20 among many others) included exchange effects (I > 0) with 
a variety of potential functions for more accurate predictions. Much of this 
work is discussed in Johnstons^® review. Available predictions for the 
species H, H 2 and He are discussed below. 

2. 3. 1. 1 Hydrogen Atoms 

The prediction of the neutral Bremsstrahlung cross-section for 
hydrogen atoms has received considerable theoretical analysis over the 
years because of its importance in astrophysics. Detailed predictions, 
including exchange effects have been presented by John^9^ Geltman20^ 

Doughty and Fraser^l and Stilley and Callaway. The results of these 
various treatments are substantially in agreement with the results generally 
being presented in tabular and/or graphical form. 

The first relatively accurate, (in comparison to detailed numerical 
calculations), analytical expression for the H atom neutral Bremsstrahlung 
cross-section was prese»>ted by Dalgarno and Lane. In this analysis the right 
side of (Eq. 2. 3) was expanded to first order in terms of the zero-order elastic 
scatteriiig cross section qQ (E), where E is the electron energy. The 
latter quantity was expanded at low energies by use of c'flective range* 
theory^^ with the result 
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(E) = 4r (Dj + D 2 E 


+Dj Ein(E) + D 4 E + DgE^^^)a 


where is the Bohr radius and the quantities may be related to molecular 
parameters. Their resulting expression for the absorption cross section 
was 

-34 -3/2 -a r 

= 4. 51 X 10 T ' (1 -e ) IDj Ij (a) 

1/2 3 . 

+ b ' D 2 I 2 (a) + b D 3 (a) + b J (a) In (b) (2. 7) 

3 3/2 

+ bD^ - Ij (a) + b ' D 5 I 4 (a) + .... ) 


where a = hV/kT, b = 1.2395/X (Wtn) and the quantities Ij (a) are integrals 
which were tabvilated vs a. These latter quantities are reproduced in 
Table 2. 1. 


Dalgarno and Lane applied Eq. (2. 7) to a ntunber of atoms and 
molecules. For the case of H atoms the relevant parameters are Dj = 11.24, 
D 2 = 7.2, Dj = 4 . 959 , D 4 = -11.58 and D 5 = -0.4023. Eq. (2.7) was shown 
to provide good agreement with the more detailed calculations of Geltman. ^0 
Dalgarno and Lane^^ suggest that for the case of H atoms relationship (2. 7) 
should be accurate for values of X and T such that (a + 2) T< 30, 000°K. 

25 

Very recently Stallcop has developed a very simple analytic 
expression for the H atom neutral Bremsstrahlung cross section valid in 
the limit where hv/kT is small. In this work the absorption cross section 
is again defined in terms of the elastic scattering cross-section, and 
after a number of mathematical approximations it is shown that 


-45 -hv/kT I 

Q 2.96 X 10 T (1 - e )"‘ (E«/hv) 

enn 

5 

• exp (- ^ (k^))/k^ , cm 


( 2 . 8 ) 


1 /2 

where k_ = (kT/E ) , Eu is the ionization potential for hydrogen atoms 

T XI 

and 
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TABLE 2.1 


PARAMETERS FOR PREDICTION OF NEUTRAL 
BREMSSTRAHLUNG CROSS SECTIONS (DALGARNO & LANE, 1966) 


a 

^l(a) 

^2 (a) 

^3 (a) 

4 (a) 

0. 02 

1. 000 X lo'^ 

1.175 X 10^ 

7. 752 X 10^ 

2. 056 X 10^ 

0.04 

2. 500 X 10^ 

2.077 xio"^ 

8. 390 X 10^ 

1.817 X 10^ 

0. 06 

1. Ill X 10^ 

7. 537 X 10^ 

2.261 X 10^ 

4. 397 X 10^ 

0.08 

6. 249 X 10^ 

3. 671 X 10^ 

8. 863 X 10^ 

1. 606 X 10^ 

0. 10 

3. 999 X 10^ 

2. 101 X 10^ 

4. 270 X lo'* 

7. 356 X 10^ 

0. 15 

1. 777 X 10^ 

7. 622 X 10^ 

1. 121 X 10^ 

1. 780 X lo"^ 

0.20 

9. 992 X 10^ 

3. 711 X 10^ 

4. 297 X 10^ 

6. 501 X 10^ 

0. 30 

4.437 X 10^ 

1. 345 X 10^ 

1.093 X 10^ 

1. 572 X 10^ 

0.40 

2.492 X 10^ 

6. 536 X 10^ 

4. 069 X 10^ 

5. 743 X 10^ 

0. 50 

1. 592 x 10^ 

3. 732 X 10^ 

1. 869 X 10^ 

2. 629 X 10^ 

0.60 

1. 104 X 10^ 

2. 359 X 10^ 

9. 805 X 10^ 

1. 388 X 10^ 

0. 70 

8. 092 

1. 599 x 10^ 

5. 637 x 10^ 

8.091 X 10^ 

0. 80 

6. 181 

1. 141 X 10^ 

3.466 X 10^ 

5. 068 X 10^ 

0. 90 

4. 871 

8.470 

2. 243 X 10^ 

3. 354 X 10^ 

1. 00 

3. 935 

6.483 

1. 511 X 10^ 

2. 319 X 10^ 

1. 50 

1.692 

2. 301 

3.091 

5.601 

2. 00 

9. 550 X lO"^ 

1.096 

9. 088 X 10"^ 

2. 047 

2. 50 

6. 032 X 10“^ 

6. 141 X lO"^ 

3. 167 X 10“^ 

9. 404 X lo"^ 
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(2.9) 


^(k^) «4.862 (1 - 0.2096 k^ +0.0170 k^^ - 

0.00968 kj^) 


valid for hv/kj < 1 , 0. 1 < k^ < 1 . 

Stallcop has shown that this^^imple relationship is in good agreement 
with more detailed predictions. Note that relationship (2. 8 ) is valid 

at considerably higher temperatures than the expression of Dalgarno and 
Lane^^, £q. (2. 7). although the latter is to be preferred at the higher values 
of hv/kT. Absorption cross sections for H atoms, as determined from 
Eq. (2, 8 ) are shown in Fig, 2. 1 vs. temperature for wavelengths of 5 and 
10 Mm. The cross sections scale as the square of wavelength and approxi- 
mately linearly with temperature over this temperature range. 


2. 3. 1.2 Hydrogen Molecules 

The neutral Bremsstrahlung absorption cross section for H 2 can be 
of considerable importance at low temperatures. For example, at a pressure 
of 10 atm the equilibrium concentration of H 2 will be greater than that of 
H at temperatures below 4100°K. Unfortunately, little theoretical effort has 
been directed towards the prediction of this cross section. 

23 

Although Dalgarno and Lanes' analysis is strictly only appropriate 
for atoms, it has been applied to the Hydrogen molecule under the assump- 
tion that the low energy elastic scattering of an electron by a molecule is 
dominated by the spherically symmetric component of the molecular 
potential field. From Ref. 23 the relevant parameters to be used with 
Eq. (2. 7) are Di = 1. 6456, D 2 = 5. 8492, D 3 = 0. 86147, D 4 = -2. 6502, 

D 5 = 0. 0. The relationship is expected to be reliable for wavelength and 
temperature such that (a + 2) T < 10,000°K. 

The only other prediction for H 2 appears to be that of Somerville, 

His analysis is similar to Dalgarno and Lands' with a slightly different 
approximation for qp (E). The latter work'^^ is preferred. However, given 
the lack of experimental data it is difficult to determine its reliability. 


2. 3. 1.3 Helium 


The most recent predictions for He have been performed by 
Geitman^^’ who used a full partial wave analysis, Eqs. (2.3) - (2.5), with 
model potentials that include a polarizability term. These predictions are 
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0 2 4 6 8 10 


T, °K X 10“^ 

Fig. 2, 1 Predicted Hydrogen Atom Neutral Bremsstrahlung 

Absorption Cross-Section vs. Temperature for X = 5, 10 
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in good agreement with earlier work by McDowell et al and are 
estimated by Geltman to be accurate to within 30%. 


Predictions of the neutral Bremsstrahlung cross sections for He 
are tabulated in Ref. 16 over the wavelength range of 0. 5 - 20 Mm and 
temperature range from 500 - 20, OOO^K. Representative values at 
X = 2. 0, 5. 0 and 10. 00 |ini are plotted vs temperature in Fig. 2. 2. In 
this wavelength/temperature range the cross section scales approximately 
with the cube of wavelength and linearly with temperature. The cross 
sections are roughly a factor of 3 lower than those predicted for H atoms. 

Theoretical predictions for a number of other atoms and molecules 
are also available. These include the rare gases, heated air 

species, 27-29 and more recently the alkali atoms^®"^2^ jhe 

work in this area is voluminous, the above references being representative 
rather than inclusive, and will not be reviewed here. 

2.3.2 Comparison Between Theory and Experiment 

There does not appear to be any experimental measurements of the 
cross-sections for infrared neutral Bremsstrahlung with H, He, H2. 

There are a limited ntimber of measurements for rare gas and air species and 
these are listed in Table 2. 2 along with theoretical predictions taken from 
Refs. 16 and 28. There is minimal overlap between the conditions of the 
various measurements. However, where there is such overlap, the measure- 
ments are in agreement with the exception of the cross section for neon. 

What is remarkable is that in the infrared, experimental measurements 
for all gases shown are consistantly higher than the theoretical predictions, 
typically by a factor of two to four. 

It is hard to rationalize this large difference between theory and 
experiment even though the experiments involved are admittedly complex. 

It would appear that considerable more effort, both experimental and 
theoretical, must be extended in this area before low temperature neutral 
Bremsstrahlung absorption cross sections may be predicted with a high 
degree of certainty. It should be emphasized, however, that the 
theoretical efforts are most appropriate to the hydrogen atom. 

2.3.3 Example - Low Temperature Seeded Hydrogen Plasma 

It is worthwhile to briefly discuss the dominant absorption process 
for a specific configuration. Equilibriiun species concentrations^^ for 
pure hydrogen at 10 atm total pressure are shown in F^g. 2. 3. Note that at 
the highest temperature shown (10,000 K) the equilibrium ion mole fraction 
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Number Density, part/cc 



3 - 





TABLE 2. 2 


COMPARISON OF EXPERIMENTAL MEASUREMENTS AND 
THEORETICAL PREDICTION FOR NEUTRAL 
BREMSSTRAHLUNG CROSS SECTIONS 


Gas 

■Sil 

■H 

IIK8RRiH!9llii 

lEHRRHRHSn 

He 

5000 

in 

• 

o 

-40 

5. 5 X 10 (a) 

-40 

3.2 X 10 

Ar 

5000 


-40 

1.8x10 (a) 

-40 

1.8x10 (g) 


9700 

3. 1 

1.0 X lO’^^(b) 
-38 

7. 5 X 10 (c) 

2.5 X 10'^®(g) 



5.0 

3.6 X 10"^^(b) 

1.0 X 10"^^(g) 

Ne 

12,600 

3.1 

30 

1.9 X 10' (b) 

6. 9 X 10 ^®(c) 

-3ft 

1. 8 X 10 (g) 



5.0 

-38 

7. 2 X 10 (b) 

-3R 

7. 2 X 10" (g) 



9.85 

-37 

7. 8 X 10 (c) 

5.2 X 10"^^(g) 

Xe 

8200 

2 

3 X ’0'^®(b) 

- 3ft 

1. 5 X lo' (g) 



3. 5 

9.6 X 10 (b) 

6.5x 10'^®(g) 



5 

2. 0 X 10“^^(b) 

1.7x lO'^^g) 
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TABLE 2,2 (cont. ) 


N 

9175 

2. 

- 38 

3. 2 X 10 (d) 

6.4 X 10”^^(e) 

•39 

8.9x10 ’(g) 

6. 9 X 10‘^ V) 



3.5 

1.3 I0“^^(d) 

2.6. 10“^^(e) 

4.0 X 10"^®(g) 

3. 1 X 10~^®(h) 



5.0 

3.0 X 10* ^^(d) 

1. 1 X 10‘^’(g) 
38 

9.8x 10"^°(h) 

O 

5000 

2.0 

2.4 X 10’^ V) 

2.7x lO'^^g) 
2.5 X lO'^V) 



5.0 

2.6 X 10‘^^(f) 

-38 

2.9x 10 ^^g) 
3. 6 X 10”^®(h) 


9700 

3. 1 

-38 

7.6 X 10 (d) 

_*o 

7.8 X 10 (c) 

_«o 

2.3x10 

1.0 X lO'^V) 


- 

5.0 

2, 6 X 10”^^(d) 

8.5x 10'^%) 

-38 

3.6x10 (h) 



9*85 

3.0 X 10"^^(c) 

1.0 X 10 ^^(g) 
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TABLE 2. 2 (cont. ) 


N2 

8350 

2.0 

-38 

3.8x10 (d) 

1.4 X 10‘^®(h) 



3.5 

1.3 X 10”^^(d) 

-38 

7.0 X 10 (h) 



5.0 

3.0 X 10"^^(d) 

2.0 X 10"^’(h) 

02 

5000 

2.0 

9.6 X 10“^V) 

3. 2 X 10"^^(h) 

1 


5.0 

1.0 X 10"^ V) 

•38 

3. 7 X 10 (h) 


a = ref. 33, b = ref. 34, c = ref, 35, d = ref. 36, e = ref. 37, 
f = ref. 38, g = ref. 1:6 h = ref. 28 
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is 5, 7 X 10 . It can readily be shown from a comparison of Eq. (Z, 2) 

to the neutral Bremsstrahlung cross sections shown in Fig. 2. 1 that neutral 
Bremsstrahlung will dominate ionic for ion mole fractions < 10“^. Thus, 
over the full temperature of Fig* 2.3 neutral Bremsstrahlung will be the 
dominant continuous absorption process. Even if a seed, such as Cs, were 
introduced a concentration of more than 2 x 10"^ in the cold gas wo\ild be 
required before ionic Bremsstrahlung would dominate. This large a 
concentration of Cs would produce the undesirable result of increasing the 
average molecular weight of the cold gas from 2 to 4.7 gms/mole. 

Thus, it would appear that even for a seeded system, neutral 
Bremsstrahlung would provide the dominant continuous absorption mechanism 
below 10,000 K with collisions with the hydrogen molecule dominating below 
4000^K. The remaining question is whether there will be any contributions to 
the total absorption coefficient from line absorption (i. e. bound-bound transi- 
tions). Very recently the infrared emission spectra of cesium contaminated 
air was measured^^ at a temperature of 3300^K under conditions where 
only 10% of the Cs was ionized. This measurement is shown in Fig. 2.4 with 
identification of the prominent Cs lines. The monochromator resolution is 
""0.2 lim over this wavelength range. The Bremsstrahlung prediction is a 
semi-empirical extrapolation from the measurements of Taylor and 
Caledonia. As can be seen there are several prominent Cs lines which 
could enhance the total absorption coefficient, although not at CO 2 laser 
wavelengths. 

40 

Hyman and vonRosenberg have demonstrated an interesting feature 
concerning the spectral positioning of these atomic lines. This is that since 
infrared lines involve pairs of levels with high principal quantum numbers 
the spectrum will be fairly hydrogenic with the exception of a relatively 
small number of lines characteristic of the atom of interest. These hydro- 
genic lines typically involve states with the same principal quantum number 
n (unallowed transitions in hydrogen) or at most involve a change in n of 
unity. Hyman and vonRosenberg point out that the line positions of these 
transitions, which will occur in most atomic gases, may be specified by 
the relationship 

2 2 

X (^m)= 0.0911 ” y— (n > 5) (2.10) 

and, thus, the first few li. .s will appear at 4. 05, 7.46, 12. 36 and 19.05 gim. 

1 he lowest of these lines would explain the ubiquitous 4 pm line which has 
been observed in a number of gases (see for instance thc‘ data in Refs, 34, 36). 
Note that prominent lines at these wavelengths are obsi rved in the Cs 
spectrum shown in Fig. 2.4. 
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RADIATIVE INTENSITY 
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Fig, 2.4 Infrared Emission Spectrum of a Cs-Air Mixture. Line 

Identification corresponds to Cs Atomic Transitions, Condi- 
tions as shown. Taken from Ref, 40. 


- 18 - 




The prediction of the absorption coefficients for such bound-bound 
transitions is beyond the scope of the present work. It shovild be stressed 
that as the temperature is raised the Cs will become completely ionized in 
which case bound-botmd transitions of the bath gas would be dominant. 

2.4 SUMMARY 

Theoretical predictions for the cross sections of inverse ionic and 
neutral Bremsstrahlung absorption have been reviewed. For the most 
part die discussion has been limited to the prediction of the infrared 
absorption cross sections for light gases, such as H, H 2 and He, at 
temperatures below 10, OOO^K. Although there is no data available for 
H, Hz and He the measured infrared cross sections for other gases are 
typically a factor of two to four higher than theoretical predictions. 


The particular case of a seeded hydrogen plasma has been discussed. 
Under equilibrium conditions the dominant continuous absorption process 
will be neutral Bremsstrahlung at temperatures below 10,000 K and seed 
levels < 2%. Hydrogen molecules will be the dominant collision partner 
b«‘low 4000K and hydrogen atoms above. At the lower temperatures 
bound-bound transitions could provide the dominant absorption process at 
specific wavelengths. The Bremsstrahlung absorption coefficient is 
linearly proportional to electron density and thus, even in the presence of 
a seed will decrease exponentially with decreasing temperature at tempera- 
tures below 2500-3000°K. 
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3. MOLECULAR ABSORPTION 


3. 1 INTRODUCTION 

As mentioned in the previous section, a basic disadvantage in using 
inverse Bremsstrahlung for laser gas heating is that under equilibrium 
conditions the Bremsstrahlung absorption coefficient decreases exponentially 
with decreasing temperature at low temperatures (T< 3000°K in the case of Cs, 
which has the lowest ionization potential of all the elements) and is negligibly 
small at room temperature. An alternate absorption mechanism which can 
be as (or more) effective in a cold gas as in a heated one is direct absorption 
by the vibration/rotation (V/R) bands of a seed molecule. This mechanism 
however does have a number of disadvantages not encountered in inverse 
Bremsstrahlung. 

1) A given molecule will have a finite number of 
absorbing transitions which will occur over a 
limited wavelength range; 

Z) The absorbed energy is deposited in the vibra- 
tional rather than translational modes of the 
molecule and efficient kinetic mechanisms are 
required to convert the absorbed energy into heat. 

In particular, when the upper state of an absorbing 
transition is populated more rapidly than it can be 
de-activated, gas "bleaching" can occur; 

3) As the gas temperature increases undesirable 
chemical phenomena, such as dissociation, can 
degrade the absorption process. 

This last point is particularly important for larger molecules. For 
example, SF^,, which is an efficient absorber of CO 2 laser radiation, will be 
significantly dissociated at temperatures above ~1700“K. Even a relatively 
stable triatomic such as CO 2 will dissociate at temperatures above ''2500°K. 
Generally speaking, only diatomic molecules with large dissociation energies 
can provide significant laser gas heating with minimal chemical effects. As 
will be discussed below, the first two disadvantages may be minimized by 
performing the absorption in a high pressure system. Because of the 
complications introduced by chemical effects, the major portion of the 
discussion below will be limited to absorption by diatomic molecules. How- 
ever, much of what is said will also hold true for larger molecules. 
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3. 2 DIATOMIC ABSORPTION 


The pressure broadened absorption coefficient for a given V/R 
transition with vibrational levels v - 1 v and rotational levels J J i. 1 
is defined by 


“v-l,J,AJ 


(V), cm 



J 


^ v-l,J,AJ \,v-l,J,AJ 
4t? c a V. 


(3.1) 




where v is wavenumber, c the speed of light, \ wavelength, A the 
Einstein coefficient of the transition, N the populations of the upper and 
lower states, g the degeneracy factor and the full linewidth at half 
height in cm“^. It can be seen from the first bracketed term on the right side 
of Eq» (3« 1) that the absorption coefficient decreases as the upper state 
population approaches the lower. This is the phenomenon of “bleaching** 
alluded to earlier. The number density dependence of the absorption co- 
efficient is of some interest. The line width is directly proportional to 
number density and thus on line center (v = "^v,J,6j) absorption co- 
efficient is independent of the total number density. On the other hand on 
the wings of the line (v AT ^ absorption coefficient scales 

as the square of the number density. Thus, there can be a large variation 
in the pressure dependence of the absorption coefficient over the line shape. 

For any given vibrational transition there are a large munber of 
allowable rotational transitions. For example, the room temperature 
absorption spectrum for the v = 0 1 transition of CO is shown in Fig. 3. 1 

It has been asuumed here that the pressure is such that the line width is 
small compared to the spacing between lines (^3.8 cm*l for CO). The 
v = 1 -> 2 absorption spectrum will have a similar shape but will be shifted 
to somewhat longer wavelengths, as shown in Fig. 3. 1, because of the 
molecular anharmonicity. As can be seen there is a considerable overlap 
between the v = 0 1 and v = 1 2 bands (and also with higher bands not 

shown) however if the line widths are small compared to the line spacing 
there will be none but accidental wavelength coincidences between the various 
lines. Thus laser light of a given wavelength can typically be absorbed by 
only one of the manifold of possible absorbing transitions. 
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4.7 

WAVELENGTH, MICRONS 


Room temperature, low pressure absorption spectra of the 
V = 0 -» 1 band of CO 


This constraint will not apply at sufficiently high pressures where 
the line width is > to the line spacing, ( ^ 30 atmospheres for self- 
broadening in CO). In this limit the molecule will absorb continuously over 
the wavelength region of the V/R band. The absorption coefficient at any 
given wavelength is no longer given directly by Eq. (3. 1), but rather now 
involves a sum of Ihe absorption coefficients of a number of overla|ping 
transitions, and this sum is independent of total number density depending 
only upon the mole fraction of the absorbing gas. \ 


An approximate band model has been developed for the prediction 
of this high pressure absorption spectra. At pressures at which the C^nes 
just overlap, occuring when = Z Be, where Be is the rotational constant 
of the molecule, it is assumed that the absorption spectra becomes co'i- 
tinuous with a magnitude at any wavelength proportional to the line center 
value of the absorption coefficient of the transition nearest that wavelength 
(from Eq. (3. 1)). The prqportionality constant is determined by the re>- 
quirement that the integration of the absorption coefficient over wavenu^'iber 
be proportional to the bandstrength, i. e. , ^ 


N r 

--wS- 


(V) dv 


\ 

(3. 2)\ 




where S is the bandstrei^th in cm" 


for an absorber number density N^, 


Nq is Loschmidts munber and N is the number density of the absorber. 


The resulting absorption spectra is defined at file discrete points v 
by the relationship 


v,j,ar 


v-l** v 


»vl-.y _ 2 

~ N . v-1, J,A J'^v, v-1, J, 


v-1 


h (2J + 1) J (J 1) he 

l6lTkT e*P - kx 



(3.3) 


Where i is the vibrational population of the absorbing state, is the 
vibrational correction to rotation, AT takes on the values 1 (A = 0 ) and 
the Einstein coefficient A^ v-l»J# A J related to that for the vibrational 
transition, A^, v-1* Hozde-London factors. It can be seen that the 

absorption cross section ^v-1, v independent of the total number 

density and varies with temperature as the rotational population. Similar 
band models have been shown to provide a good approximation^^ to the actual 
high pressure absorption spectra. 
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This band model has been applied to the molecule CO and predictions 
for this case are shown in Figs. 3. 2 - 3.4. In Fig. 3.2 the room 
temperature absorption cross section vs wavelength is shown for the lowest 
three vibrational transitions of CO. The cross section increases with in- 
creasing vibrational level because of the variation in the Einstein coefficient 
(a linear dipole moment function for CO was emfdoyed in the calculations). ^3 
The important feature is that, because of the large overlap between the 
absorption spectra for each vibrational transition, laser light of a given 
wavelength may be absorbed simultaneously by a number of vibrational 
levels. This has the effect of minimizing bleaching phenomena; for now, 
if a particular transition approaches the bleaching limit the upper state 
of that transition can also absorb the incident laser radiation. 

Of more interest from the viewpoint of laser heating is the 
temperature variation of the absorption cross sections. This behavior 
will, of course, be quite sensitive to the choice of incident laser wavelength. 
Predictions for the absorption cross sections of the first five vibrational 
levels of CO over the temperature range of 200 - 300U^K are shown in 
Fig. 3.3. The incident laser radiation is taken to be at a wavelength of 
4. 755 t-im, which would correspond to a P branch v = 1 0, J = 10 CO 

laser. Although a high power CO laser of this type has not yet been 
developed, the results may be taken as generic for any diatomic molecule 
undergoing absorption at a wavelength near the center of its fuiuiamental V/R 
band. As can be seen, above 300 K the absorption cross sections for the 
various vibrational levels are all of the same order of magnitude. The 
absorption cross sections for individual transitions generally tend to de- 
crease at the higher temperatures shown because of the effect of the 
rotational partition function. However, absorption on the higher vibrational 
transitions (v = 4 5 as well as higher transitions not shown) peak in the 

higher temperature region because absorption occurs in the far rotational 
qiiantum niunber wings of these bands. The net effect of this is that the 
overall molecular absorption coefficient can be relatively insensitive to 
temperature variations depending, of course, upon the vibrational distri- 
bution function. 

A similar res^llt is shown in Fig. 3.4 for an incident laser wave- 
length of 4.855 ^m. This wavelength corresponds to an attainable P branch 
CO laser transition, v = 3 2, J = 8. Since this wavelength occurs in the 

wing of the CO v = 1 band, see Fig. 3.2, the absorption cross section 

for this transition initially increases rapidly with temperature. This is 
an important point since the upper vibrational states of CO are not 
significantly populated at room temperature. However, it is clear that as 
soon as ground state molecules are excited to the first vibrational state 
through absorption, excitation to higher states (including transitions to 
states higher than those shown in Fig. 3.4) will be quite rapid. 
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4.6 4.8 5.0 

wavelength, Mm 


Fig. 3.2 Room temperature, high pressure absorption spectra of the 
V = 0-»-, V = l-»2andV = 2-»3 bands of CO. 


V, v+l 



Fig. 3. 3 Absorption cross section vs temperatures for the five lowest 
vibrational transitions in CO at a wavelength of 4. 755 ^m. 
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It would appear then that optimum absorption by diatomic 
molecules can be realized in the high pressure regime characterized 
by overlapping lines, since in this limit laser light of a given wavelength 
may be absorbed simultaneously on a number of different vibrational 
transitions, thus minimizing bleaching. This multiline absorption should 
also have the effect of flattening out the temperature variation of the total 
molecular absorption coefficient. It now remains to discuss techniques for 
the rapid conversion of the absorbed vibrational energy into translation. 

3. 2. 1 Conversion of Vibrational Energy to Translation 

The molecule CO was chosen as an example not only because its 
molecular properties are well defined, but also because it has a very high 
dissociation energy, 11.1 eV, and thus may be heated to relatively high 
temperatures (~ 6000°K) without introducing significant chemical effects. 
One disadvantage of this molecule however is that the rate constants for 
vibrational deactivation are quite small at temperatures below 1500°K 
(see for example Ref. 44). These small rate constants tend to preclude 
efficient coupling between the vibrational and translational modes. There 
are however, several atomic species such as Fe, H and O which have been 
shown to be efficient (k = 10"^^ - 10"^^ cc/sec) vibrational deactivators of 
CO;*^^"^^ however, it is difficult to introduce such species int«.) a cold gas 
flow. 


One technique which has been considered which, although not a 
cyclic system, may be of some interest for systems concerned with 
development of thrust, etc, , involves the use of a rare gas seeded with 
iron pentacarbonyl, Fe (00)5. Fe (CO^ ) absorbs in the same wavelength 
region as CO, its vibrational modes will strongly couple with translation 
and, when dissociated, produces both the absorbing molecule CO and the 
efficient deactivator Fe. Unfortunately, detailed modeling of this system 
is complicated by the fact that the low temperature absorption spectra of 
Fe (00)5 is not well defined. 

To properly model gas laser heating by diatomic molecular 
absorption one reouircs a computer code which includes the phenomena of 
multiline absorption, vibration to translational (V-T) deactivation and 
vibration-vibration ( V- V) cxchajige between ail vibrational levels. The 
resulting vibrational distributions will not only be cut of equilibrium with 
translation but can also be non-Maxwellian. 

Preliminary calculations have been performed for CO neglecting 
V-V exchange processes and assuming that the absorption occurs on only 
one -rotational transition in each vibrational band. In this instanc*', the 
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master equation for the time variation in vibrational population density of 
a constant density gas may be written 


d Ny/d t 


I 

hv 


v-1 -» 


(V) 


1 ” 


- a T A, N 
V, J, AT 


■1 


I 

hv 


V~* V + 1 




|nv- 


+ 1 , 


J', AJ* N ^ 

V + 




k. , M. ^ N - exp ( ^ i i 

i,v-»v-l,v i|v kT / v-1 J 

c. , , M. \ N , - exp ( ^ ^v + I, vl ( 

1 , v+l-v x\ v+1 i^--/ \ ] 


where 1 is the laser intensity, v i?= laser frequency, AE^ v 1 
difference in vibrational energy be. een levels v, v-1 and ^ is the 

rate constant for species Mj de-activating level v to v-1. The quantities 
ay^ are defined as 


%,J,-1 ~ (~ P branch (3.5) 

and 

/ 2 (j + 1) he B~\ „ , , /, /. 

a T . 1 “ ( ~ » R branch (3.6) 

v,J,+l \ kT / 

where J is the rotational number of the transition closest to the laser wave- 
length. 


A calculation was performed for conditions corresponding to an 
initial pressure of 30 atm with 0.2% CO and 0.04% Fe in a helium bath. 

The laser was assumed to have an intensity of 10^ watts/cm^ at a wave- 
length of 4.855 U m, corresponding to the absorption cross sections of 
J'ig. 3.4. The time histories of both gas absorption length and temperature 
for this case are shown in Fig. 3, 5. As can be seen, the gas temperature 
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Gas temperature and absorption length vs time. Conditions 
as shown. 
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ABSORPTION LENGTH. CM 


30 ATM 99. 75% He/0. 2%CO/0. 04% Fe 



increases essentially linearly from 400 to 3000^K and^ more important, 
the absorption length remains approximately constant over this temperature 
range* The initial drop in the absorption length is caused by the early time 
increase in the populations of the upper vibrational states. The inclusion of 
V-V exchange in the modeling is not expected to produce a major change in 
these predictions. 

3.2.2 Other Molecules 

There are, of course, a number of candidate diatomic molecules 
which can be de-activated relatively efficiently by stable species. For 
example, the room temperature rate constant for self-vibrational deactivation 
of NO, i*e., M. NO, is cc/sec (Ref. 48, among others) and 

furthermore the absorption sp< '^'.ra for NO is quite similar to that for CO. 

The predicted absorption spectra for the lowest two vibrational transitions 
of NO are shown in Fig. 3 .6. The middle peaks in these spectra arise 
because the NO ground state is a Tf state and thus the Q branch is allowed. 

(For simplicity in this calculation a linear dipole moment was assumed and 
the small wavelength shift between spin states was neglected). Shown at the 
top of the figure are the wavelength positions of several readily attainable CO 
laser lines. It would appear that NO would prove to be a viable absorber for 
CO laser radiation. 

It should be emphasized that this technique may be generally applied. 

In general, the radiation from most high power gas lasers which operate on 
vibrational transitions can be absorbed by several different diatomic 
molectiles if the total gas density is sufficiently high to produce overlapping 
lines. The one major exception to this statement appears to be the 10. 6^m 
CO 2 laser. Although a number of the diatomic metal oxides can absorb 
radiation of this wavelength, these species are not readily produced under 
the conditions of interest. Of course this same line broadening phenomenon 
occurs for larger molecules* These could be reviewed to determine which 
species could optimally absorb CO 2 laser radiation wdth minimal chemical 
effects. 

3. 3 absorption; by larger molecules 

There will of course be a large number of molecules, triatomic or 
larger, capable of absorbing laser light of any specified wavelength. Only 
two examples of these will be discussed below. The first of these will be the 
relatively large molecule SF^,, which has been employed in many systems 
as an absorber of CO 2 laser radiation. The second will be the molecule 
CO 2 , which is taken as an example of a relatively stable triatomic. 
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The species SF£, is an example of a large molecule whose vibra- 
tional kinetics and infrared absorption properties have been examined in 
some detail* SF6 has six fundamental vibrational modes which, including 
degeneracies, correspond to fifteen vibrational states. Vibrational excita- 
tion in any mode is rapidly redistributed among the manifold of vibrational 
states, the room temperature constant for this process being ~ 3 x 10* ^ 
cc/sec for collisions with rare gas atoms. Vibration to translation de- 
activation then occurs predominantly from the lowest energy vibrational 
mode at a somewhat lower rate, the room temperature V-T rate constant 
for collisions with the rare gases being approximately 0. 6 x 10"^^ cc/sec.^^ 


The absorption properties of SFf, at 10. 6 kJim have received con- 
siderable study at temperatures near 300°K (for example. Ref. 50 and 
references therein). is somewhat unusual in that the vibrational mode 

which absorbs the 10. 6 (Jim radiation consists of an estimated 1000 lines/ 
cm*l and at pressures above 0. 5 torr the absorption becomes essentially 
continuous and independent of total pressure. The room temperature 
absorption cross section^O for SF^ is -10-17 

cm^, a factor of 100 larger 

than the peak value for CO, Fig. 3.2, and thus considerably less SF6 
seeding would be required to provide any specific amount of absorption. 


50 

Anderson has performed laser saturation measurements in SF^ 
at room temperature and 1 atmosphere total pressure and found that sig- 
nificant bleaching begins to occur at laser intensities of 3 x 10^ w/cm^. 
(This would correspond to saturation laser intensities of 10^ w/cm^ at 
pressures of ^^30 atm.) Presumably the saturation intensity woiald increase 
with increasing temperature since it has been shown^^ that the rate constant 
for vibrational deactivation of SF^ increases as . 


The temperature dependence of the SF6 absorption coefficient is not 
well known and is difficult to evaluate because of the increased importance 
of ’’hot band** absorption, i. e. , absorption by excited states, with increas- 
ing temperature. Anderson^O has measureu the absorption coefficient over 
the temperature range of 300 - 360° K and found that it decreased linearly 
by approximately 20% over that range. A linear extrapolation of his data 
would imply that the absorption coefficient goes to zero at 573°K, which is 
clearly a lower bound. However, as the temperature exceeds lOOO^K 
chemical decomposition will become important causing a gradual deteriora- 
tion of the absorption process. 


Thus at room temperature SF(^ absorption at 10. 6 m would appear 
to be as effective as the CO/Fe system at 4.855 m which was discussed 
earlier and has the additional advantage that the number of absorbing 
molecules which must be introduced into the flow is considerably lower. 
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Unfortunately, gas heating to temperatures much larger than 1700°K is pre- 
cluded because of thermal dissociation of the SF^. 

The vibrational kinetic behavior for SF^ discussed above is common 
for many larger molecules. Energy absorbed in any given vibrational mode 
is first distributed among all the available vibrational modes and then con- 
verted to translation, primarily through vibrational deactivation of the lowest 
energy vibrational mode. For example, this behavior is observed in CO 2 
although the rate constants^^ for these respective processes are lower than 
in SF^. 


The two fundamental vibrational bands of COz are centered about 4. 3 
and 15 tJim, The 15^ band is a relatively weak absorber. However, the band 
at 4. 3 |Jim can absorb more strongly than CO. The same type of high 
pressure absorption analysis that was applied earlier can be used to predict 
the temperature dependence of the CO 2 absorption coefficient; however, the 
required band model would be considerably more complicated. Absorption 
by these bands would deteriorate at temperatures > 2500 K because of CO 2 
dissociation. However, this temperature is sufficiently high to cause thermal 
ionization of a seed species such as Cs. 

CO 2 can also absorb at 10. 6 in a "hot band" transition, the in- 
verse of the CO 2 laser transition, however, the absorption coefficient for 
this process can vary significantly at low temperatures. Recently, ^4 this 
absorption coefficient was measured over the temperature range of 300 - 
1600°K, for the P (18) and P (20) laser lines, at pressures where the line 
widths were considerably less than the line spacing. The absorption co- 
efficient increased from a value of 2 x 10“^ cm"^ at 300°K to a peak value 
of 3 X 10"^ cm“l at 800°K anH then decreased monotonicallv to 10"^ cm"^ 
at 1600°K. (All of these absorption coefficients are for an absorbing molecule 
number density of Nq. ) The temperature variation between 300 - 800'^K is 
determined primarily by the vibrational population of the lower state of the 
absorbing transition and thus would be similar even if the pressure were 
such that the lines overlapped. 

There are obviously a manifold of such molecules which can be 
considered. However, the general phenomenology of interest from the 
viewpoint of laser gas heating should be adequately represented by the two 
examples discussed above. 

3. 4 SUMMARY 

The conversion of laser energy to gas translational energy via 
absorption by the vibration/ rotation bands of molecules has been examined. 

The most effective mechanism would appear to be absorption by the V/R 
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bands of diatomic molecules at sufficiently high pressures so that the 
lines of neighboring transitions overlap. Beth the effects of gas bleaching 
and the temperature variation of the absorption coefficient may be 
minimized in this high pressure region. Unfortunately, there does not 
appear to be any readily produceable diatomic molecule which can absorb 
CO 2 laser radiation. There are a number of larger molecules which can 
absorb 10. 6 ^m radiation; however, efficient high temperature absorption 
by these species is precluded because of thermal dissociation. Such 
species may be sufficient, however, to raise the gas temperatures high 
enough for inverse Bremsstrahlung to be effective (~2500^K for a Cs seed). 

In the kinetic analysis presented above, the bath species were 
assumed to be inert gases such as helium. If H 2 were the bath gas, addi- 
tioxial kinetic mechanisms, both vibrational and chemical, would have to 
be considered. 
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4. PARTICULATE ABSORPTION 


4. 1 INTRODUCTION 

Gas laser heating may also be achieved through seeding the gas flow 
with particulates. In this case, the particles absorb the laser radiation, are 
heated to temperatures above that of the bath and then thermally conduct 
heat to the surrounding gas. The advantages of this technique are three- 
fold; (1) the absorption can be spectrally continuous, allowing use of any 
laser wavelength, (2) the absorption process is efficient in a cold gas and 
only a small percentage of the absorbed energy remains in the particle, 
and (3) the conversion of the absorbed energy to gas translational energy is 
relatively rapid. The basic disadvantages are (1) the more complicated 
system required for reliable particle formation and injection, and (2) 
particle vaporization at increased temperatures. 

In the sections below, the phenomenology of particle heating will be 
described. The discussion includes a review of the absorption properties of 
particulates and preliminary modeling of the gas heating rate. 

4. 2 DISCUSSION 

4. 2. 1 Particulate Absorption 

Large particles, defined as particles of characteristic size much 
larger than the wavelength of the incident radiation, will absorb as gray 
bodies, with an absorption cross section defined as 

cJa (^) - 

where is the wavelength dependent material emissivity and R is particle 
radius. For simplicity, the particles are assumed to be spherical. As the 
particle size approaches the wavelength of the incident radiation, the 
particle causes a perturbation to the electromagnetic field about it, and 
the full Mie theory^^ is required to predict the absorption cross section. 

In this instance 0^ complicated function of particle radius and the 

complex index of refraction, M = n - in* of the material. A simplification 
in Mie theory results when the particle size is considerably smaller than 
the wavelength of incident radiation (Rayleigh absorption). Under these con- 
ditions the absorption cross section is given by 
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This type of absorption is referred to as volume absorption and generally 
written as 


°A = “ '" • 


(4.3) 


where V is the particle volume and a is the volume absorption coefficient 
defined by 


a 



- 1 
+ 2 


( 4 . 4 ) 


The maximum absorption efficiency per unit mass loading of 
particulates is achieved in the limit of volume absorption since in this 
limit the total volume, rather than surface area, of the particle is involved 
in the absorption process. It is clear that volume absorption will apply if 
the particle is optically ”thin**. This condition is satisfied when 


a R < 1 (4. 5) 

The value of a can vary considerably depending on material, v/ave- 
length and temperature. Generally speaking, the volume absorption 
coefficient for elemental materials such as C, Al, etc. will have a 
relatively smooth spectral variation whereas the absorption by more 
complex materials such as Si02» AI2O3, etc. can be highly structured 
spectrally. An example of the infrared spectral variation of a at room 
temperature for two such materials is shown in Fig. 4 . 1 . The data for 
soot is taken from the work of Volz^^ while that for Si02 is taken from 
Steyer et al. In the visible the values of a for carbon are approximately 
a factor of 6 higher^® while on the other hand Si02 is essentially transparent. 

The characteristic absorption length for particles in the volume 
absorption limit is given by 


L = (4/3 ttR^ Npa) 


( 4 . 6a) 
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or 


L = p /(M a) (4.6b) 

H H 

where Np, Mp and are particle number density, mass loading and 
specific density respectively. Thus, to achieve an absorption length of 
10 cm for nominal values of a = 2 x 10^ cm*^, Pp = 2 gms/cc would 
require a particle mass loading, Mp, of 10“^ gms/cc. Furthermore, from 
condition (4. 5), the particle radius would have to be less than 5 ^ m. 

This nnass loading would correspond to 1% of the typical reservoir densities 
under consideration and to a particle density of Np 10”5 /r^ (for example, 
10*^ particles/cc for a radius of 1 ^m). 

4. 3 HEATING PHENOMENOLOGY 


Upon absorption of laser radiation, the particle will heat up, 
thermally conduct heat to the surrounding gas, melt and vaporize (the small 
effect of thernaal radiation is neglected). The effect of vaporization is 
only important at higher temperatures, approaching the boiling point of 
the material, and will be neglected in the following analysis. 


For the case of a constant density gas, the differential equations 
defining the temperature history of a laser irradiated particle and the gas 
surrounding it may be written 


4/3 TT la = 4/3 TT p„ d (Cn T)dt 

F Fp p 


2 1 


(4.7) 

- 4 TT R kd To/dr 

P IR 

r = R 


— 1 



T )J / dt -V (k VT) = 0 

r > R 

(4.8) 


where vaporization and melting have been neglected and the temperature field 
within the droplet has been taken to be spatially uniform. In these equations, 
I is the laser intensity, Cpp is the particle specific heat, k is the gas 
^ermal conductivity, Cv is the gas specific heat at constant volume and 
T is the average gas temperature. If the gas is seeded with Np particles/cc 
of radius R, the rate of increase of average gas temperature is given by 


pd (Cy T)/dt = - 4n R^ Np k d T/ dr 



(4.9) 
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Unfortunately, the quantities a, Cp , and k are themselves 
functions of temperature and it is not possible to derive analytic solutions 
for Eqs. (4. 7)-(4. 9). Nonetheless, one can examine the qualitative 
behavior of these equations by examining their solution when the former 
quantities are constant. 

Under this assumption, one can readily show that the solution to 
Eq. (4, 8) for r *** R is given by59 


whence 


T - T — (T - T) , 
r p 


dTp/Srj =-(T„-T)/R 


(4.10) 


(4.11) 


If relationship (4. 11) is substitu^"' ’ into Eq. (4.7), it can be shown 
that the solutions of Eq, (4, 7), (4. 9) are 


Tp - T = (Tp - T) 


j^l - exp (-t/T ) J 


(4. 12) 


and 


dT/dt = 


L-' I 


p Cv (1 +FCo /C 


Pp' 




1 - exp (-t/T ) 


] 


(4. 13) 


where F is the relative mass loading, Mp/p, and 


(T - T) = ^ ^ ■ ' 

' P ^s.s. 3 k (1 + F Cpp/Cy) 


(4. 14) 


T = 


^P ^Pp ^ 


3k (1 + F C„ /C ) 
Pp V 


(4.15) 


In essence the difference between the particle and average gas temperature 
approaches a steady state value, Eq. (4. 14), in a characteristic irradiation 
time T, Furthermore, the gas heating rate approaches a constant value 
with the same time constant T, After this time, the absorbed laser energy 
is converted almost totally into gas heating. 
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There are a number of constants which must be specified in 
evaluating Eqs. (4. 12)-(4. 15). Room temperature values for several 
representative absorbing materials and bath gases are given in Table 4. 1. 
(Cpp, and k all increase with increasing temperature. The largest 
variation is in Cpp which in the case of carbon increases by a factor of 
"“3 between 300 and 2000 K). It can be seen from these representative 
values that the quantity F Cpp/C^, which is found in Eqs. (4, 13)-(4. 15), 
will be significantly smaller than unity for the relative mass loadings of 
interest, 10*^, and can be neglected. 

The efficiency of gas heating _by particulate absorption may be 
optimised by minimizing both (Tp - T)s^s. can be seen from 

Eqs. (4. 14)“(4. 15) that this may be accomplished by reducing the particle 
size. For a specific numerical example, if I = 10^ w/cm^, a = 2 x 10^ cm-1, 
and k, Pp, Cp took on characteristic values then a particle of size 
0. 3 i-Lm would^approach steady state in a characteristic time of 5 ^ sec 
with a temperature differential, Tp - T, of 500^K. These latter two 
quantities would decrease by an order of magnitude if the particle size 
were reduced to 0* 1 |Jim. 

It is particularly important to minimize the temperature differential 
since the above relationships break down as the particle temperature 
approaches the boiling point. In this limit the major p>ortion of the 
absorbed energy is lost in the latent heat of vaporization rather than in 
gas heating (a similar, but transient, effect occurs at the melting point). 

More importantly, as the particle vaporizes and thus reduces in size, 
the absorption coefficient decreases. Finally, if the bath gas is not inert, 
undesirable heterogeneous chemical reactions may occur at the particle 
surface. The boiling temperatures for several absorbing materials are 
listed in Table 4. 1. Vaporization effects, however, may be important at 
temperatures considerably below these values. 

Interestingly enough, although the particle absorption coefficient 
will decrease as vaporization occurs, the vapor species may provide 
alternate absorption mechanisms . For example, at temperatures > 2500°K, 
the dominant metal oxide vapors will be diatomic. Absorption in the 
V/R bands of these molecules occurs at ^ 10 m, i. e. , in the wavelength 
region of the CO 2 laser. Furthermore, the vapor species, or their elemental 
components, can have ionization potentials significantly lower than that of 
the bath gas (see Table 4. 1), thus allowing equilibrium inverse 
Bremsstrahlung absorption to dominate at lower temperatures than in the 
pure bath gas. 
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TYPICAL PARTICULATE PROPERTIES 
( T • 300® K ) 
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To summarize, although detailed predictions require a full 
computer solution, it has been demonstrated that particulate absorption 
can provide a viable method for high efficiency laser gas heating. The 
technique is optimized by utilizing very small particles, < 0.3 iim 
with typical mass loadings of order 10"^ gms/cc. Although the effects 
of vaporization are deleterious to particulate absorption, the vapor species 
themselves can provide alternate absorption mechanisms and it is 
conceivable that a system design which would allow "complete" vaporiza- 
tion of the particles could be very efficient. The basic disadvantage of 
particulate seeding is the equipment requirements for particle production. 
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5. ANALYSIS OF FLOW AND STABILITY REQUIREMENTS 


5. 1 INTRODUCTION 

This section presents the part of our program that was directed 
toward fluid dynamic considerations for laser heated flows. Specifically, 
the objective was to develop the capability to predict the thermodynamic 
and fluid dynamic properties resulting from the radiation heating of a 
nozzle flow* In addition, the capability to determine the stability of the 
resulting flow field is considered essential to a study of the efficiency 
of converting laser energy to a high velocity flow* 

For simplicity, a quasi one -dimensional nozzle flow is considered 
(Fig. 5. 1). In this nozzle geometry the laser radiation will enter either 
through a window at the flow inlet plane or through the exhaust at the exit 
plane* The flow is considered to be that of a perfect gas, and the 
absorbed laser energy will be considered to be instantaneously converted 
to temperature (i. e. equilibrium flow)* The equations of motion will be 
solved to obtain steady state solutions, and a **local'* stability analysis will 
be performed to determine how stable the absorption zone is to small 
disturbances. 

The stability of radiatively heated flows along with suggestions for 
acoustic wave amplification by radiative absorption have been examined to 
some extent by many authors. In general, the technique for deter- 
mining stability is to obtain a steady state solution and then super-impose 
an infinitesimal disturbance. The disturbance is composed of a number 
of discrete partial fluctuations, each of which consists of a wave. For 
each single oscillation of the disturbance the sign of roots to the resulting 
characteristic equations are examined as a function of intensity and wave 
niunber to determine when stable solutions can he expected. Since the 
steady flow properties are not constant in a one-dimensional nozzle flow, 
when waves are introduced, the corresponding boundary value problem 
is rather difficult. However, much can be learned about wave propaga- 
tion without solving the full boundary value problem, and information can 
be extracted from a study of the modified dispersion relation. To obtain 
the modified dispersion relation in inhomogeneous media the actual 
disturbance is approximated by harn^ionic functions over a small distance 
(small compared with the characteristic lengths associated with ibe 
gradients of the flow) and, hence, the analysis is only an accurate 
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description "locally"* This procedure is commonly used in problems of 
this type (e. g. Monsler^^ has examined the stability of radiative heating 
between two infinite plates). 

5.2 STEADY NOZZLE FLOWS WITH LASER-ENERGY ADDITION 

Consider a one-dimensional flow through a nozzle of variable area 
ratio that is heated by laser radiation (Fig. 5. 1). The laser beam enters 
either through the nozzle inlet plane in the same direction as the flow 
(termed Option A) or through the exhaust plane in the direction opposite 
to the flow direction (termed Option B). The governing equations for 
quasi-one -dimensional steady flow without viscous dissipation, diffusion 
and thermal conduction, but including absorption of laser radiation are: 



where p is the density, u is the velocity. A is the area, p is the pressure, 
h is the enthalpy, I the local laser intensity and is the absorption co- 
efficient. The minus and plus signs correspond to Option A and B, 
respectively. 

Calculations will be carried out for generalized absorption coeffi- 
cients of the form 

k = K p" T*" 

V o 

where Kq is a constant, T is the gas temperature, and n and m are 
constants characterizing the absorber. It is useful to write the governing 
equations in terms of an optical depth T defined as^2 
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k d X 
V 


(5.2) 


T 


X 


Inserting Eq, (5,2) into the radiative transfer equationflast of Eq. (5, 1), we 
can immediately integrate the resulting equation to obtain the laser power 
distribution (I A) as a function of optical depth T, The result for Option A 
is 

- T 

I A = Ij A^ e 


where the subscript i denotes conditions at the nozzle inlet plane, and 
for Option B 


T-T 

IA = I^A^e 


where the subscript e denotes conditions at the exhaust exit plane. When 
T is introduced into the remaining equations, and the dependent variables 
are made dimensionless with 0j» u., hg..,,etc. , where i denotes inlet 
conditions, Eqs, (5,1) become 



(5.3) 


lE 

d T 


+ p u 


d u 
dT 


0 


(5.4) 



(Option A) 


(Option B) 


(5, 5) 


The important parameter F, which is the ratio of laser power to the initial 
total energy flux of the flow, is defined as 
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where IP is the incident laser power (equal to liAj for Option A and 
le Ae for Option B) and hg^ is the initial stagnation enthalpy. Since 
Eq, (5. 5) is directly integrable, the increase in stagnation enthalpy, 
hs = h + u^/2, becomes 


hg = l+ r^l-e ^ (Option A) 

-'’’e / T V 

hg = l+ re ^e -1/ (Option B) 


(5. 6) 


It can be seen from Eqs, (5. 2) and (5. 6) that without proper absorption co- 
efficient tailoring, the absorption will take place close to the inlet or exit 
plane. The actual physical location will depend upon the chosen nozzle 
configuration and can only be determined by inverting Eq. » 5. 2). 

The equation of state 


P = PRT 


together with Eqs. (5.3) - (5.4) provide the equations for the remaining 
variables P, u. and T. A general method for solutions to this system of 
equations may be employed. The governing equations can be combined 
to yield a relationship which governs the Mach number variation through 
the nozzle 


dM^ 

dT 


-F(M ,T) 

l-M^ 


(5. 7) 


and 





- , dh 

'hT df 


2 ±A 
A dT 


where M is the Mach number. For convenience, the nozzle contour is 


-4Q. 



chosen to be the one parameter nozzle configuration used by Buonadonna 
et al, 64 




(I -t-d -T) e^- (1 +d) 
- (1 f d) 


(5. 8) 


where Aj is the initial area (T = 0) and A»p is the throat area located at 
T = d* The calculations are not restricted to this configuration and can 
easily be extended to any arbitrary nozzle shapes* As a matter of fact, 
the nozzle contour will be one of the parameters to be examined when 
optimizing the power conversion efficiency. 

Eq. (5. 7)together with Eq. (5.6) and (5.8) are used to obtain the 
Mach number distribution along the nozzle axis. Since we are only in- 
terested in a flow which accelerates from subsonic (convergent section) 
to supersonic (divergent section), the solution is obtained by firs.^* locating 
the position corresponding to the singular point when F (M^, T) 0 as 
1, and then integrating upstream and downstream to get the desired 
Mach ntimber distribution. 

Example calculations for area ratio A'p/ H = ?. and 
r= Ofl, 2 and 4 are obtained for two nozzle shapes as shown in Figs. 5.2 
and 5. 3. It can be seen that when heating is in the supersonic flow 
region (Option B) the flow is driven toward M = 1. This is because the 
nozzle area ratio does not increase rapidly enough in the supersonic 
region to compensate the laser heating effect. It can be shown that a faster 
increase in area ratio in the divergent section provides a continuous in- 
crease in Mach number distribution. The above results are presented in 
terms of the optical depth T, However, they can easily be transformed 
to the actual physical space for given absorbers by inverting Eq. 2). 

Knowing the Mach number profile for a nozzle configuration with 
known inlet conditions and laser power, the complete steady state solu- 
tion for p, u, and T, can be calculated. Further examples of such 
calctilations are discussed in Section 5.4. 

5. 3 STABILITY ANALYSIS 

Due to the inhomogeneous nature of the one-dimensional nozzle 
flow, it is difficult to solve the resulting eigenvalue problem when 
perturbations are superimposed upon the steady solution. Hence, the 
”locaF^ stability analysis described in Section 5. 1 will be applied to deter- 
mine the stability at each point along the nozzle. 
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Fig. 5. 2 Mach number distributio 
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Since the superimposed wave disturbance will be time dependent, 
it is necessary to write the full unsteady one -dimensional equations with 
heat addition. They are 


£p b (pu) Pu d A _ 

&t & X A d X 


(5.9) 


* du , „ 5 u , ^ - 


(5. 10) 




dp _ ^ ^ ^ 

p ~ P T 


(5. 12) 


r — = k = K p T 

ox V o 


(5.13) 


Only Option A has been considered; however, the analysis can be easily 
extended to Option B. For simplicity, the pressure is eliminatec^ from 
Eq. (5. 10) and (5. 11) by means of Eq. (5. 12), Now we have 4 equations 
for 4 unknowns (p, u, T and T), and can proceed to superimpose disturbances 
upon the previously obtained steady flow field by assuming 


P = pQ (x) + P' (X, t) 


u = u^ (x) + u' (x, t) 
T = (x) + T* (X, t) 


T = T (x) + T' (x,t) 


(5. 14) 


where P^, T^, Uq, and are the known steady state variables and the 
primes indicate small perturbations, i. e, , p’ <Kp^ Neglecting the second 
or higher order terms and substituting Eq, (5. 14) into Eqs. (5.9) - (5. 13), 
yield a linear set of equations for the perturbation quantities. 
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(5.15) 
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It can be seen in Eqs« (5* 15) that the coefficients are only functions of the 
steady state solutions and hence, for the sake of simplicity, the above 
equations are non-dimensionalized with respect to the steady state variables, 

i. e. , 


P 



t 

u 



o 


X 



and t 


t u 
o 

L 


where the characteristic length L is taken to be 1/k j . With this normaliza- 

*^o 

tion, Eq, (5, 15) becomes: 
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where Y is the ratio of specific heats, a = R and 


r = r 



-T 

o 


e 


The steady state flow field is assumed to be influenced at a given location 
by a disturbance which is composed of a number of discrete partial 
fluctuations (harmonics). The perturbations representing a single 
oscillation of the disturbance are assumed to be of the form 
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(5. 17) 


F-'?(n. K)e' <K»-nt) 


— '^/o K'l ^ -fit ) 

u = U (u, K) e 


T = x'ta K) 

T ='t‘ (n,K) 


The nondime ns ional wave number K is a real quantity and X = 2 TTL/K is 
the wavelength of the disturbance* The quantity, Q, is the complex 
dimensionless frequency 


n =: flj. + 

where 0^. is the frequency of the partial oscillation and 

(the amplification factor) determines the degree of amplification or 
damping* The disturbances are damped if fli< 0 and the steady -state 
flowfield is stable, whereas for Qi > 0 disturbances will grow. The velocity 
of propagation of the wave in the x-direction (phase velocity) is given by 


C = 



u 

o 


Substituting Eq. (5. 17) into Eqs. (5. 16), we obtain 


(~Jtnu + A - in + iK^ + (-^^nPo + a - An A + 1 k) 

'dx ° dx ' 

( An « + a An T + ia k)^ + {-^ An Uq - ifl +i K 
Ux ° dx ° 'dx / 

+ (a An p + i aK) t' = 0 

' d X ° ' 


(5.18) 


+ Z An u - ifl + i k) u + (- 1’* m - i —.0 + i a k) 1^+ = 0 

Sx o ' \ y-1 Y-1 / 


- 1^- ^ 


iK T = 0 
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The above eqiiation can be written as 


A V = 0 


where A is a 4 x 4 matrix and V is the unknown column vector* i. e. * 


V 



-1 


In order that the solution exists* the determinant of A must vanish. Hence 


-^n «o -It fn A-in +iK ^ fn + — £nA + iK o 

iSx Sx dx o g- 


dx 

u„ + a Jtn T + iaK 

^ dx o 

-~Xnu_ -iQ+i K 
dx ° 

a-“AnP_ + kxK 
dx ® 

0 

Y-1 

An T + — An u 

dx ® Sx 

a -^=-An T +2-^*n\i_ - f* m 
Y-1 dx o dx ^ 

r' = 

- r' n 

+ to.0 

-in+iK 

K 


- n 

0 

-m 

iK 


(5.19) 


This expression is the modified dispersion relation and provides the relation- 
ship between the real Kand the complex Cl • The method of solution to 
assume K and* then* solve for the complex 0 (3 roots) at each location along 
the nozzle. This corresponds to applying an infinitesimal wave of wavelength 
2TTL/K at each location and examining whether the disturbance will grow* 

(Qi > 0) or decay (0^< 0). The locus of = 0 define contours of neutral 
stability. The states of neutral stability can be one of two kinds. For 
Or = 0 the amplitudes of a small disturbance can grow or be damped 
aperiodically (stationary mode) when Clr 4 0 (oscillatory mode) the theory 
tells us not only what the mode of disturbance is that will be manifested at 
onset of instability* but also the characteristic frequency of the oscillations. 
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5.4 RESULTS AND DISCUSSION 


Sample solutions to Eq. (5. 19) have been obtained for Option A (laser 
beam from upstream) and for one nozzle configuration (d = 2 and A/A^ = 0. 5). 
These calculations can also be carried out for Option B as well as for other 
nozzle shapes. 

For this particular nozzle configuration calciilations have been done 
for r= 0, 1, 2 and 5. As shown in Fig. 5.4, the Mach number distributions 
are plotted as functions of , the steady- state flowfield optical depth. For 
a given absorber (and hence a known absorption coefficient) one can readily 
transform these results into the jdiysical space by inverting Eq. (5.2). In Figs. 
5. 5 - 5. 7 the density, the velocity and the temperature distributions versus 
are presented. As can be expected, the acceleration of the flow near the inlet 
region, where the absorption takes place, increases with increasing laser 
power. As the laser power is increased, the temperature also increases and 
the density decreases such that the pressure profile becomes a weak function 
of r. For these steady- state flowfields disturbances were superimposed and 
solutions to Eq. (5. 19) were examined to deternnine whethe. such disturbances 
grow or decay in the nozzle. 

Stability maps for F= I, 2 and 5 are presented in Figs. 5. 8 - 5. 10. 
These results are for the generalized absorption coefficient of the form 


k^ p T (5.20) 

and a sample case where the absorber is fully ionized such that n = 2 and m = -3/2 
is computed. Therefore, absorption is via inverse Bremsstrahlung. For this 
case the gas must first be heated to a temperature high enough to produce a 
sufficient number of electrons for this absorption process to dominate (unless 
a discharge is initiated via gas breakdown). 

Since the governing equations have been normalized with respect to the 
steady- state solution, the above results are only fxinctions of F (ratio of laser 
power to the initial steady energy flux), the nozzle shape and the generalized 
form of absorption coefficient. 

In Figs. 5. 8 - 5. 10 a nondimensional wavenumber is plotted against the 
optical depth in the steady state flowfield, and contours of (1^ = 0 are drawn. 

These contours form the boundary between stable and unstable regions in wave- 
number space. At each nozzle location (corresponding to a specific Tg) the 
shaded zone indicates the range of wave numbers which result in growing dis- 
turbances. For example, in Fig. 5. 8 at = 1 disturbances with wavenumbers in 
the range 0. 1 ^K< . 37 and 6 < K_< 4. 5x 10 will grow; however, at = 1, 7 dis- 
turbances with wavenumbers over the entire range will be damped. It must be re- 
membered that this analysis will not tell us whether a disturbance, which is initiated 
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Fig* 5*4 Mach number distribution along the nozzle 








Fig. 5, 5 Dimensionless density 
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Fig. 5. 6 Dimensionless velocity distribution 
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Dimensionless temperature distribution 
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Fig. 5, 8 Contours of neutral stability for T = 1 
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Fig. 5. 9 Contours of neutral stability for f = 2 
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OPTICAL DEPTH. \ 


Fig, 5, 10 Curves of neutral stability for F = 5 
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and grows in one region, will dampen after it travels to a new location. 

The analysis serves as an indicator as to where potential absorption wave 
phenomena may be initiated. It is desirable to design a nozzle/absorber 
arrangement which contains a minimum number and extent of these unstable 
zones. 


The unstable zones are numbered 1, 2 or 3, corresponding to the 
three roots of 0. Often at a specific location only one or two of these roots 
yields an unstable zone. For example, in P*ig. 5. 8 at = . 3 the second 
root results in the unstable region for 0.1^ K < .4 and the first and third 
roots yield stable waves. At = 1 the first and third roots result in tin- 
stable zones, but the second root yields stable waves. In addition, one root 
may yield several zones of instability. For example, in Fig. 5. 9 the third 
root yields two distinct unstable zones, one at.2<To<.5, .9<K<1. 1 
and the other at . 8 < Tq < 1. 3, . 1 < K < , 8. The unstable regions are 
upstream of the throat and near the inlet area as expected where the 
absorption takes place. It can be shown that the regions labeled 2 and 3 
are physically unimportant by inverting Eq. (5. 2) and returning to physical 
dimensions. This can be seen if we refer to Fig, 5, 11. 

Figure 5, 11 presents the stability map translated into physical 
coordinates, x is the distance along the nozzle axis and most of the laser 
radiation is absorbed from x = 0tox = 10 cm. For simplicity, the sample 
gas has been chosen to be helium with 0. 1% of cesium to provide a fully 
ionized seed that absorbs by inverse Bremsstrahlung at the initial 
temperature Tj = 5000® K. The initial pressure pi is 5 atm and as stated 
the cesium concentration Csj. = 0. 1%, The unstable zones are shaded and 
are located from x = 0tox = 10 cm in the region where the absorption takes 
place. It is now necessary to determine which anticipated disturbance wave 
numbers are physically important. We must first determine an appropriate 
length scale to compare with the disturbance wavelength. The steady state 
solution for this example case indicates that by x = 10 cm most of the laser 
power has been absorbed. Therefore, a characteristic absorption length 
scale is 10 cm. If a disturbance haj a wavelength much greater than 10 cm, 
the absorption will only affect a fraction of the wave and growth of the 
entire wave due to unstable heating will not be significant. Therefore, only 
disturbances of wavelengths less than 10 cm can be considered important. 
Since K = 2 TT/X, where X is the disturbance wave length, this implies that 
only wave number K >2 TT/10 cm"^ = . 6 cm“^ are of interest. Unstable 
region #1 is the only important region. It is now desirable to make 
this region as small as possible to minimize the possibility of unstable 
heating in the nozzle. As the laser power is increased this unstable 
region decreases in size as is seen in Figs. 5.9 and 5.10. Therefore, a 
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Fig. 5. 11 Curves of neutral stability 
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nozzle design based upon the conditions in Fig. 5. 10 may be considered a 
good design for stable conversion of laser energy to kinetic energy. 


5. 5 SUMMARY 

Analyses of the steady flow resulting from laser heating in a 
nozzle configuration and the resulting stability of such flows have been 
presented. It must be emphasized that the sample results of Figs. 5. 1 - 
5. 11 are not necessarily for the most efficient nozzle design to convert 
laser power to high velocity exhaust. The analyses must be applied to 
other nozzle shapes and absorbers in order to optimize Ihe power con- 
version efficiency. 
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